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Abstract 


This report first provides some background information on intrinsic point 
defects, and on carbon and oxygen in silicon in so far as it may be relevant 
for the efficiency of solar cells fabricated from EFG ribbon material. We 
discuss the co-precipitation of carbon and oxygen and especially of carbon and 
silicon self interstitials. A simple model for the electrical activity of 
carbon-self-interstitial agglomerates is presented. We assume that the 
self-interstitial content of these agglomerates determines their electrical 
activity and that both compressive stresses (high self- interstitial content) 
and tensile stresses (low self-interstitial content) give rise to electrical 
activity of the agglomerates. The self-interstitial content of these 
carbon-related agglomerates may be reduced by an appropriate high- temperature 
treatment and enhanced by a supersaturation of self- interstitials generated 
during formation of the p-n junction of solar cells. 

It is suggested that oxygen present in supersaturation in carbon-rich 
silicon may be induced to form SiC> 2 precipitates by self-interstitials 
generated during phophorus diffusion. It is proposed that the SiO_-Si 
interface of the precipates gives rise to a continuum of donor states and that 
these interface states are responsible for at least part of the 
light-enhancement effects observed in oxygen containing EFG silicon after 
phosphorus diffusion. 


•Permanent address: Max-Planck-Institut fur Metallforschung, D-7000 Stuttgart 

80, West Germany. 



1. Introduction 


Silicon ribbon material grown by the edge-defined film- fed growth (EFd) 
technique invariably contains a high concentration of carbon, about 10^cm“3. 
The carbon is introduced by the slotted graphite die in which liquid silicon 
rises by capillary action. Silicon solidifying above the die 

is withdrawn in the form of a thin continuous ribbon. EFG ribbons 
usually exhibit a much smaller minority carrier diffusion length L D (typically 
in the order of 10-50 pm) and also a correspondingly lower minority carrier 
life time than float-zone grown (FZ) silicon or Czochralski-grown (CZ) 
silicon. Recently it has been suggested that the presence of carbon-point 
defect agglomerates may be at least partially responsible for low values in 
EFG ribbons. 1 Other and more commonly involved causes for low values in EFG 
ribbons are the presence of a comparatively high density of crystallographic 
defects, such as various types of grain boundaries and dislocations (both in a 
very inhomogeneous distribution), as well as contamination with metallic 
impurities. Grain boundaries in EFG ribbons can be excluded as the main cause 
for low L D values, because the average distance between electrically active 
grain boundaries is usually in the order of hundreds of pm and therefore much 
exceeds L D . Therefore, we are left with a "diffusion-length triangle", see 
Fig. 1, with the remaining three possible causes for the low L D values in EFG 
ribbons (dislocations, metallic impurities, and point-defect-carbon 
agglomerates) denoted at each corner. These three types of defects may 
influence each other (e.g., dislocations may absorb point defects; 
point-defect agglomerates and dislocations may act as gettering centers for 
metallic impurities; metallic impurity precipitates may act as nucleation 
sites for point-defect agglomerates; etc) and in different areas of the 



ribbons different types of defects may be responsible for a low L Q value. 

The influence of metallic impurities and of dislocations on L D has been 
treated in other publications. We therefore concentrate in the present report 
on the possibility of a diffusion-length limitation by point-defect-carbon 
agglomerates. The report is by no means meant to be complete or to provide 
final answers. It is rather our intention to stimulate more specific 
experimental investigations by proposing some speculations and also a very 
rough model for understanding at least qualitatively the change of L D by 
various tretments of the ribbons, such as high temperature heat treatments of 
phosphorus diffusion steps. Since our model is based on interactions between 
intrinsic point defects and carbon and oxygen in silicon we will first compile 
some basic information on intrinsic point defects (Section 2) and on carbon 
and oxygen in silicon (Section 3) which is required to understand and to 
develop our model in some detail in the subsequent sections. Readers who are 
familiar with point defects and carbon and oxygen in silicon may directly 
proceed to Section 4. 


2. Some basic information on point defects in silicon. 


2.1. Point defects and diffusion under thermal equilibrium conditions. 

A perfect crystal can decrease it Gibbs free energy by introducing a 
certain concentration of point defects into the crystal. For a crystal 
consisting of one element, like silicon, two basic types of point defects have 
to be considered, namely vacancies, V, (one silicon atom missing) and 



self-interstitials, I, (one extra silicon -atom squeezed into the crystal) as 
schematically shown in Fig. 2. In Fig. 2, we have also indicated chemical 
point defects, such as substitutional dopants and foreign interstitial atoms, 
such as most metallics in silicon. The thermal equilibrium concentrations (in 
dimensionless atomic fractions) of vacancies, C®^ and of self interstitials, 
C® q , are given by 

= exp (sj/k) exp (-H*/kT) (1) 

= exp (sj/k) exp (-Hj/kT) (2) 

where k is Boltmann's constant and T the absolute temperature. and are 
the entropy and the enthalpy of formation of a vacancy, respectively. s| and 
are the corresponding quantities for self-interstitials. 

Although in the literature various values can be found for H^, these 
quantities have not yet been determined experimentally, basically because the 
concentrations of vacancies or self-interstitials are very low (< lO^cnT^) 
even at the melting point. In the literature there exists a long standing 
controversy whether self- interstitials are present at all in silicon in 
quantities comparable to vacancies. 2 In the following we take the view that 
vacancies and self-interstitials co-exist in silicon^ and have to be taken 
into account for understanding diffusion processes in silicon. (Gosele and 
Tan^, Frank et al. 1 *). Whereas hardly anything reliable is known on C® q and 
the situation is different for the basic silicon material transport 
quantity, namely the self-diffusion coefficient D^, which is given by 



4 


where and Dj are the diffusion coefficients of vacancies and 
self-interstitials, respectively. The individual terms in (3) are given by 
3,4 

=0.57 exp (-4.03ev/?) cm 2 s" 1 
D I C i q = 914 exp (-4; 84 eV/?) cm 2 s _1 , 
as also shown in Fig. 3. 

The individual diffusivities and Dj are as poorly known as the 
concentrations and C® q . We, nevertheless, present in Fig. 4 and Fig. 5 

recent rough estimates of these quantities by Gosele and Tan^, but remind the 
reader that this information has been reached very indirectly and may have to 
be revised considerably in the future. We mention explicitly that the 
quantities expressed in (1-5) hold for intrinsic silicon and depend on the 
position of the Fermi level. 

2.2 Interaction between vacancies and self-interstitials 

If both vacancies and self-interstitials are present in silicon it can be 
expected that they react with each other according to the reaction 


(4) 

(5) 


I + V = 0, 


( 6 ) 



where 0 represents the undisturbed silicon lattice. This reaction, which for 
silicon has first been proposed by Prussin^ and later on more specifically by 
Sirtl?, involves the recombination of I and V as well as the spontaneous 
thermal creation of I and V pairs (Frenkel pairs) in the silicon lattice. If 
the point defect concentrations are disturbed by some external means (e.g. by 
surface oxidation, which injects interstitials into the silicon) then for 
sufficiently long times reaction (6) will establish a local dynamical 
equilibrium between I and V which may be described by the law of mass action 


c i c v * c ? qc ? q 

In eq. (7) Cj and Cy are the actual concentrations of self- interstitials and 
vacancies, and and C| q their respective thermal equilibrium values. 
Measurements of Antoniadis and Moskowitz® showed that at 1100° C it takes a 
time 


T dyn (H00°C) * 1 hr (8) 

to reach local dynamical equilibrium conditions described by (7) . It has been 
estimated^ that T dy>n should increase with decreasing temperatures 
approximately inversely proportional to the silicon self-diffusion coefficient 
D sd . For 900°C this leads to the estimate: 

T dyn (900 ° C) = days ' (9) 


which means in practical terms, that for experiments performed at 900°C at 
much shorter times (e.g. a 30 min phosphorus diffusion) the vacancies and 



self-interstitials can be treated as non-interacting species. 


\ 2.3 Diffusion of substitutional elements 


The diffusion of substitutional elements in silicon (such as Group-Ill or 
Group-V dopants, substitutional carbon, or as a special case, silicon atoms 
themselves) requires intrinsic point defects. Under thermal equilibrium 
conditions, their diffusivity D s is composed of a contribution D y involving 
vacancies and a contribution D® involving self-interstitials as "diffusion 
vehicles" according to3»^»8»9 

D s = + Dj (10) 

The relative vacancy and self-interstitial contributions depend on the 
specific element and on temperature. Generally, atoms smaller than silicon, 
such as boron or phosphorus show an essential relative contribution D|/D s of 
the diffusion component involving self-interstitials. It has been argued^ 
that carbon which is very much smaller than silicon should also diffuse 
predominantly via a diffusion process involving self- interstitials; (see also 
Section 3.2) . 

If the point defect concentrations are disturbed by external means (e.g. 
oxidation, P-diffusion) and attain values C y and Cj then the diffusivity D s of 
a substitutional element is changed to 

D !x f Dl (C l/cf> * D5<C V/ C V «<1) 


( 11 ) 



At high temperatures and sufficiently long times and Cy are related via 
> (7) . At lower temperatures and not too long times and Cy may be regarded 
as independent of each other, as explained in the previous subsection. 


2.4 Generation of non-equilibrium concentrations of intrinsic point defects 

There exist various possibilities to change point defect concentrations 
from their thermal equilibrium values C® q and C® q to some other values and 
Cy, respectively. In this subsection we do not consider particle irradiation 
in which the same numbers of vacancies and self-interstitials are generated, 
but rather surface treatments such as oxidation or the in-diffusion of 
phosphorus from the surface. In the subsequent subsection we will also 
mention non-equilibrium situations related to the crystal growth process. 

Surface oxidation of silicon leads to the formation and growth of an 
amorphous SiO., film on top of the silicon. The Si0 2 film growth leads to the 
injection of silicon self-interstitials 3>8>9 into the silicon crystal giving 
rise to c j>c| q . The resulting supersaturation of self- interstitials, 

3 I = (C I " C® q )/C® q (12) 

depends on surface orientation, and on the detailed composition of the 
atmosphere and on temperature. For (100) surface orientation in dry oxygen 
one obtains approximately ^ 



(13) 


Sj = 6.6 x 10“ 9 t 1/1 ‘ exp (2.53 eV/kDa 1 / 4 , 

which means that Sj strongly decreases with increasing temperatures. For 
v (111) surfaces the right-hand side has to be multiplied by a factor of about 

O. 6. For wet oxidation multiplication factors larger than unity have to be 
applied.^ For high temperatures and long times (e.g. , for t£. 1 hour at 
1100°C) the oxidation induced self- interstitial supersaturation (Cj > c® q ) is 
related to a vacancy undersaturation (C y < c® q ) according to eq. (7). 
Oxidation-induced self-interstitials may agglomerate and form 
interstitial-type stacking faults, so called .Oxidation induced stacking 
£aults, or in short, OSF. Oxidation- induced self-interstitials also leads 
also lead to enhanced diffusion (see eq. 11) of substitutional elements with 
a sufficiently high diffusion component involving self- interstitials (e.g., B, 

P, Al, Ga, In, and to a smaller extent, As) and to retarded diffusion of 
substitutional elements which prefer to migrate via a vacancy mechanism (e.g., 
Sb) 12 . 

The addition of a clorine-containing compound in the oxidizing atmosphere 
(e.g., HC1, Clg) may lead to the injection of vacancies into the crystalline 
silicon 1 ^* 11 *. An example for this effect in terms of Sj i3 shown in Fig. 6 
as a function of the HC1 content of the atmosphere at 1200°C. A negative 

value of s^ ^ self- interstitial under saturation) corresponds to a vacancy 
supersaturation. 

The diffusion of phosphorus into silicon starting form a high surface 
concentration creates a supersaturation of intrinsic point defects which 
diffuse into the silicon. 2 ® This supersaturation gives rise to enhanced 
diffusion of substitutional dopants in front of the in-diffusing phosphorus 



leading to the well-known "emitter-push effect" and to the broadening of 
so-called marked layers, see Fig. 7. This supersaturation also leads to 
the climb of dislocations. Screw dislocation will turn into helices; an 
example is shown in Fig. 8. Analysis of the climb direction of such 
helices 2 ^ and also of dislocation loops 22 showed that phosphorus generates a 
supersaturation of self-interstitials. This supersaturation strongly 
decreases with increasing temperature. 2 ^ At 900°C C^/Cj^ can be well above 
100. If there are no sinks for self- interstitials in the silicon bulk then 
the phosphorus-generated self-interstitials may diffuse very deep into the 
silicon, e.g., for 30 minutes (of phosphorus diffusion) at 900°C, (a typical 
situation in solar cell processing) the self-interstitial may diffuse about 

Oil 

100 pm into the silicon. n The actual range depends on the concentration of 
sinks, (e.g. dislocations) in the silicon. A high dislocation density will 
limit the self-interstitial supersaturation to a region of a few microns near 
the surface. The effect of dislocations on the range of the supersaturation 
will be discussed in more detail in subsection 2.6. 

If helical dislocations are found in phosphorus diffused silicon the 
point-defect 3upersaturation may conveniently be estimated from the number of 
turns per length, N^, via the formula 2 "*: 

ln(C I /c® ( l) * 2flN L GbVkT, (14) 

where b is the length of the Burgos vector and G the shear modules. T is the 
absolute temperature at which the helical dislocations were formed. 


2.5 Agglomeration of intrinsic point defects during crystal growth 



During the growth of silicon single cry-stals the crystal is cooled down 
from the melting temperature to room temperature. Since the point-defect 
thermal equilibrium concentrations are highest at the melting point the 
' crystal has to get rid of part of this point defect when it cools down. If 
the rate v of crystal growth is sufficiently small the point defects may 
attain their equilibrium values for the lower temperatures by the diffusion of 
excess point defects to the surface or back to the melt. For higher growth 
rates the point defects will form point-defect agglomerates. For silicon 
basically two types of microdefects which are commonly termed swirl defects, 
have been detected by Cu decoration or preferential etching: A-swirls and 

B-swirls. ° A-swirls have been identified by high-voltage electron microscopy 
as interstitial-type dislocation loops. 2 ?»28 B-swirls, if not decorated by 
copper, cannot be seen or identified in an electron microscope. 
Experimentally it has been found that carbon facilitates the nucleation and 
extends the range of stability of B-swirl defects. 2 ^ >29 ,30 Although it has 
never conclusively proven experimentally, we adopt the model of Foil et al., 
29,30, in W hi C h it is assumed that B-swirls consist of co-precipitates of 
carbon and silicon self-interstitials. These carbon-self-interstitial 
complexes are thought to be able to transform into A-swirl defects. The 
alternative view that B-swirls are carbon- vacancy complexes is unlikely, at 
least based on simple volume considerations. 

With increasing growth rates first the A-swirls and then the B-swirls 
disappear. The growth rate value v where the B-swirls disappear increases 
with increasing carbon concentration. 2 ^^ 0 Increasing the growth rate even 
further, a new type of micro-defects occurs, (see Fig. 10), so-called 
D-swirls, the nature of which has not yet been identified. 3 1 it is likely that 
D-swirls consist of vacancy agglomerates.^ 2 


In recent years it turned out that it Is not the growth rate v itself, 
which determines the appearance and disappearance of the various swirl-types 
but rather the quantity 32,33*34 

q = V/G, (15) 

where G = |AT/Azl is the temperature gradient. 

Regarding micro-defects in silicon EFG silicon ribbon material the first 
impression is that none of the concepts for microdefects formation developed 
for CZ and FZ-single crystals can be applied, basically because the growth 
rate of ribbons is much higher, the ribbon thickness is very small and the 
ribbon material contains a high density of dislocations. However, a closer 
look shows that the situation concerning microdefects may not be very 
different in EFG ribbon material. Since both v and the temperature gradient G 
are higher by about a factor of ten for EFG ribbons compared to CZ-grown Si 
crystals, the physically relevant parameter q = v/G remains in the same order 
of magnitude. The fact that the thickness of ribbons (about 200-300 jjm) is 
in the same order of magnitude as the thicknress of the swirl-free zone 2 ^“33 
near the rim of CZ or FZ-grown Si crystals does not mean that no swirl-like 
defects could be expected in silicon ribbons at all. The thickness d of the 
swirl-free zone at the crystal surfaces should be proportional to the cooling 
time, or in terms of v and G: 

d 3 <x 1/ (v G) (16) 

Since v and G are both about a factor of ten larger during ribbon growth, we 
expect about 100 times smaller denuded zones, for ribbons, in the order of 



some pm at most. A similar consideration also applies to the influence of 
dislocations. As long as the density of dislocations which can act as sinks 
for point defects is so that their average distance is 10 pm or more 
v- (corresponding to dislocation densities < lO^cnf 2 ) dislocations are not 
likely to prevent swirl formation. 

Finally we note that EFG silicon ribbons contain a high concentration of 
carbon which favors the formation of B-swirl defects 2 ^ >30. Taking into 
account all the various factors concerning crystal growth of EFG ribbons we 
consider it as likely that these ribbons contain small B-swirls which may be 
viewed as carbon-rich zones decorated with silicon self-interstitials. The 
experimental result that no carbon agglomerates or precipitates could be 
detected with transmission electron microscopy^S i s in line with this 
assumption, since undecorated B-swirls in FZ or CZ crystals could also not be 
detected by transmission electron microscopy 2 ^ » 3° . We suggest to do Cu 
decoration experiments in order to show the existence of B-swirls in EFG 
ribbons. 

In Section 4 we will present a model for the expected behavior of such 
carbon-self-interstitial agglomerates under the influence of high temperature 
heat treatments or phosphorus diffusion steps. 


2.6 Absorption of point defects by dislocations 


Intrinsic point defects present in supersaturation may be absorbed by 
dislocations. If the absorption is limited by the diffusion of point defects 



to dislocations, then the average time t until the supersaturated point 
defects are absorbed by dislocations of density o is given by 

' fc a “ 1 /(PDi,V> (17) 

where or Dy is the diffusivity of the respective diffusing point defects. 
Eq. (17) may be expressed in a different form as 

TdT = 1/p (18) 

which basically says that the average distance a diffusing point defect can 
cover before it is absorbed by dislocations is about the average spacing 
between dislocations. This means, e.g. , that a density of dislocations (which 
are assumed to act as perfect sinks for point defects) of about lO^cm -2 will 
limit a phosphorus-induced supersaturation of self- interstitial to the first 
10 |jm below the surface. 

In the context of dislocations in silicon EFG ribbons the observation of 
Kastner and Hesse 2 ^ may be of interest, who found that deformation-induced 
dislocations do not absorb point defects efficiently. A heat treatment at 
1200°C for 1 hour transformed the dislocations into a configuration in which 
they acted as almost perfect sinks for point-defects. This behavior is most 
likely related to the fact that dislocations in Si directly after deformation 
are usually dissociated and get at least partly constricted only after a high 
temperature treatment. 37 Since most dislocations in EFG ribbons are probably 
not grown-in from the melt but are introduced during cooling due to high 
mechanical stresses it is conceivable that part of these deformation- induced 



dislocations do not act as efficient point-defect absorbers before an 
additional high temperature heat treatment. 


3. Some basic information on carbon and oxygen in silicon 

Information on carbon in silicon and its impact on devices is compiled in 
a recent article by Kolbesen and Muhlbauer.38 The behavior and the properties 
of oxygen in silicon have been dealt with in the review articles by PatelSS^O 
and by Gosele and Tan. 1 ^ For detailed references we refer to these papers. 


3.1 Configurations, solid solibility, electrical inactivity, and influence on 
lattice parameter 

3.1.1 Carbon 

Carbon is mainly dissolved on substitutional sites in silicon (Fig. 1 , 
Fig. 11). It is smaller than silicon and leads to a decrease of the lattice 
constant (see Fig. 12) corresponding to a lattice contraction of about one 
atomic volume for each carbon atom incorporated.^ Substitutional carbon is 
electrically inactive . The solid solubility of carbon in silicon has been 
measured by Bean and Newman 1 ^ a s 

Cfq : i| x 10 21 * exp (-2.3ev/kT)cm”3 f (ig) 

u 

see also Fig. 13. If both carbon and oxygen are present in silicon in high 
concentrations the concept of a solid solubility becomes rather ill-defined 



because of interactions between carbon and oxygen (see for example: Zulehner 


3.1.2 Oxygen 

Oxygen is incorporated in silicon in a slightly off-centered interstitial 
position between two neighboring silicon atoms ^-48 ^ indicated in Fig. 14. 
Oxygen atoms in this position (denoted as interstitial oxygen 0 i ) is 
electrically inactive. Its solid solubility in silicon is given by ^ 

C** = 1.53 x 10 21 exp (-1.03 ev/kT)cm" 3 ; (20) 

see also Fig. 13. Because of its incorporation in interstitial position, 
oxygen, contrary to carbon, leads to an increase of the lattice constant (see 
Fig. 12) . 50 


3.2 Diffusion of carbon and oxygen in silicon 

Due to its interstitial position oxygen does not require intrinsic point 

defects for diffusion through the silicon lattice. A diffusive jump of 0 i 

requires the breaking of a Si-0 bond and simultaneously leads to a change of 

the preferential Si-O-Si axis from one <11 1> direction to another one (Fig. 

14). The diffusivity D^ of interstitial oxygen by this process is given by 
48,51 ,52 

= 0.23 exp (-2.56 eV/kT)cm 2 sec" 1 


t 


( 21 ) 



see Fig. 15. 


Due to its substitutional position the diffusion of carbon in silicon requires 
intrinsic point defects. As discussed in Section 2.3 on the diffusion of 
substitutional elements in silicon the diffusivity of substitutional carbon53, 

D = 1.9 exp (-3.1 eV/kT)cm 2 sec -1 , (22) 

c 

should be composed of a contribution involving vacancies and a contribution 
involving self- interstitials: 

D c = D I + D v (23) 

Nothing is known experimentally on the ratio D^/D^ but following the arguments 
of Gosele and Tan^, we expect that carbon as a very small substitutional 
element should show a large component D^ compared to Dy (D^»D^). A probable 
diffusion mechanism of substitutional carbon (C s ) via self- interstitials (I) 
involves the formation of a highly mobile carbon-self-interstitial complex 
(Cl) according to 

C 3 + I=(CI) (24) 

Such a highly mobile carbon-self- interstitial complex has for example been 
found by Watkins and Brower^ after low temperature electron irradiation of 
silicon. The configuration of this complex is shown in Fig. 16. 



3*3 Influence of intrinsic point defects on -carbon and oxygen diffusion in 


silicon 


Following the arguments on the dominant diffusion mode of carbon 
involving self-interstitials, we expect that a supersaturation of 
self- interstitials (C>C®Q) should enhance the thermal equilibrium diffusivity 
D c of carbon to a value 

D c Xt = D c (C I /C i q) * (25) 

As discussed in subsection 2.4 a supersaturation of self-interstitials can be 
generated by thermal oxidation or by the in-diffusion of phosphorus starting 
from a high surface concentration. In the case of phosphorus diffusion, say 
at 900°C, an enhancement of the carbon diffusivity of up to a factor of 100 
can be expected. In the case of a high carbon concentration it has to be 
taken into account that the phosphorus supplies only a limited amount of 
self-interstitials and that it may take a finite time until the excess 
self- interstitials necessary to form the highly mobile (Cl) complexes are 
formed, 

Stavola et al.^2 and Newman et al.^5 have shown experimentally that a 
supersaturation of self-interstitials may enhance oxygen diffusion at low 
temperatures (< 400°C), probably by the formation of a highly mobile (0 i i) 
complex according to-’® 

°i + I = (Oj.1) , 


(26) 



The main difference to the case of carbon is that the diffusion process 
according to (26) is a parallel reaction to the normal interstitial oxygen 
diffusion process and requires very much higher values of C j/Cj q than in the 
case of substitutional carbon. We therefore expect, that contrary to the case 
of carbon at high temperatures (say>800°C) oxidation or phosphorus- induced 
self-interstitials should not appreciably enhance the diffusion of oxygen in 
silicon. 


3.l| Volume considerations concerning oxygen or carbon 
agglomeration/precipitation 

If present in supersaturation (C-xifQ) and if the temperature is 
sufficiently high for them to diffuse, oxygen interstitials will precipitate 
in the form of amorphous or crystalline Si0 2 (or an oxide SiO x with x close to 
2). Since the volume Q qx per Si0 2 is about twice the atomic volume of silicon 
in the silicon lattice a corresponding volume increase is associated with Si0 2 
precipitation. Therefore, precipitate growth can proceed only if the silicon 
matrix is plastically deformed (punching out of dislocation loops) or if for 
every two oxygen interstitials incorporated into the precipitate about one 
silicon self-interstitial is injected into the surrounding silicon matrix*. 
The supersaturation of self- interstitials resulting from the latter 

•In principle, the absorption of a vacancy from the silicon matrix is 
equivalent, provided sufficient vacancies are available in the material. In 
the following, we will restrict ourselves to a discussion in terms of silicon 


self-interstitials . 



process gives rise to the formation and growth of interstitial-type 
dislocation loops and extrinsic stacking faults. 1 * 0 >*** >56-58 j n pig. 17 the 
basic idea for these volume considerations are depicted schematically and 
compared to the case of carbon agglomeration, where just the opposite volume 
changes occur. 

If carbon is present in a concentration above its solubility limit ( 1 9 ) a 
thermodynamic driving force will exist for a precipitation process which 
should lead to the formation of silicon carbide (SiC) provided the interface 
and stress energies (or an insufficient diffusivity at low temperatures) do 
not prevent 1 the precipitation. As an alternative to SiC formation a simple 
agglomeration process of carbon without carbide formation should also be 
considered. Both processes lead to a local volume contraction of about one 
atomic volume for each incorporated carbon atom. Therefore growth of the SiC 
precipitate or the carbon agglomerate can proceed only if for every carbon 
atom incorporated a self-interstitial is absorbed form the silicon matrix 
(Fig. 16 ). This is just the opposite behavior to that of the precipitation 
of oxygen in silicon. 

3.5 Co-precipitation of carbon and self-interstitials — 

In a typical cooling situation after crystal growth silicon 
self-interstitials are present in supersaturation (CjM^) and experience a 
thermodynamic driving force to precipitate. Since in silicon the critical 
radius for the growth of interstitial-type dislocation loops is fairly high 2 ® 
the co-precipitation with carbon in three-dimensional carbon-self-interstitial 
agglomerates (supposedly the "B-swirls") is a more favorable alternative, 
provided the carbon concentration is sufficiently high. The volume 



considerations described in the previous subsection would indicate a 1:1 ratio 
of carbon atoms and self-interstitials in the agglomerate. More sophisticated 
arguments led Foil et al.^'^O as3U mption that at high temperatures 
carbon and self-interstitials (in an extended configuration 1 **^) may form a 
separate phase which may contain a higher percentage of self-interstitials in 
the agglomerate and which is stable only in the presence of a sufficiently 
high supersaturation of self-interstitials. Without this supersaturation the 
self-interstitials will have a tendency to dissociate from the agglomerate and 
to leave behind a carbon-rich zone. 

Within the model of Foil et al. 2 3»30 B-swirls may collapse and transform 
in inter3titial-type dislocation loops (A-swirls). A similar process may be 
triggered by a phosphorus- induced supersaturation of self-interstitials. Abe 
et al.^ 0 found a high density of stacking faults below phosphorus-diffused 
areas in carbon-rich material (Fig. 19). This result indicated that carbon 
or carbon self-interstitial agglomerates may possibly act as nuclei for 
stacking faults. 


3.6 Carbon agglomerates as nuclei for Si0 2 precipitates. 

It i3 well established that oxygen in high supersaturation in silicon may 
form Si0 2 nuclei and precipitate out independent of the presence of carbon in 
the silicon. For low oxygen supersaturatlon carbon has been shown to 
facilitate nucleation of Si0 2 precipitates. 61-66 This behavior is also 
expected from the simple volume considerations of subsection 3.4: SiO 

2 

formation leads to an expansion, whereas C agglomeration leads to a 
contraction of the lattice. 



It has experimentally been shown by de Kock^ that a phosphorus-induced 
self-interstitial supersaturation leads to Si0 2 precipitation in regions, 
where no Si0 2 precipitation occurred without this supersaturation; see 
schematic Fig. 19. We assume that the mechanism of facilitating Si0 2 
formation is an indirect one and involves a sufficiently high carbon 
concentration*: The P-induced self-interstitials in supersaturation enhance 
the diffusivity of substitutional carbon according to eq. (25) which in turn 
facilitates the formation of small carbon-agglomerates, which then will act as 
nuclei for Si0 2 precipitates. For more details see ref. 68. 


3.7 Electrical activity of Si0 2 precipitates 

At temperatures around 450°C, oxygen in silicon forms small agglomerates 
which act as shallow donors ("thermal donors ") ^ . Here we are not 

concerned with these "thermal donors" which may easily be annealed out at 
temperatures above about 600°C, but rather with the electrical activity of 
real Si0 2 precipitates. It has been found that Si0 2 precipitates- formed 
between 600° and 900°C are also electrically active and associated with donor 
states. ^ *69-72 These precipitate-related donors, which show a continuous 
distribution of states, are commonly termed "new donors" in order to 
distinguish them from the lower temperature "thermal donors" mentioned above. 
It is assumed that these donor states are related to the Si0 2 -Si interface of 

* The carbon concentration of the Si sample used in the experiments has 
unfortunately not been measured. 



the precipitate and probably result from incomplete oxidation (excess 
silicon) in the SiO^ precipitate near the interface. These donor interface 
states may be destroyed by a high temperature heat treatment which allows the 
interface to relax and to get rid of excess silicon. 

The experimental results on the influence of carbon on these "new donors" 
may be summarized as follows^ 1 >69,70. carbon itself is not responsible for 
the electrical activity, but it is helpful in providing nucleation sites for 
the Si0 2 precipitates, the interfaces of which are the real sources of these 
donor states. 

The electrical behavior of these "new donors" is so similar to the 
centers found in oxygen containing EFG ribbons^ j 74 after P-diffu3ion and 
showing the so-called "light-enhancement effect"* that it is very tempting to 
asume that these oxygen related light enhancement centers are actually nothing 
else than the "new donors" found in CZ-silicon (see also Sect. 5). 


4. A qualitative model for the influence of carbon on the minority carrier 
diffusion length in EFG ribbons. 

•This means basically that with increasing photon flux the shallow 
donor states get progressively occupied, so that the electrical activity 
decreases and the minority carrier diffusion length increases with increasing 
photon flux (light intensity). 



4.1 Basic assumptions 


The essentials of our model are contained in the following four 

assumptions: 

i) As grown carbon-rich EFG silicon ribbons contain a high density >10^cm“3) 
of carbon-self- interstitial agglomerates. 

ii) For a given number of carbon atoms in such an agglomerate its electrical 
activity depends on the content of self-interstitials in the agglomerate. 
The electrical activity has its minimum (or may be not present at all) for 
a self-interstitial content corresponding to a basically stress-free 
state. A higher as well as a lower self-interstitial content leads to 
stresses as the agglomerate-matrix interface and to a correspondingly 
higher electrical activity (Fig. 20) . . 

iii) The self- interstitial content of such agglomerates may be changed by 

high-temperature heat treatments or by a supersaturation of 

self-interstitials induced by phosphorus diffusion starting from a high 
surface concentration. As a reason for a change in the self- interstitial 
content rather than in the carbon content we invoke the much higher 
diffusivity of silicon self-interstitials (Fig. 4). 

iv) Increases (decreases) in the electrical activity of these 
carbon-self-interstitial agglomerates result in correspondingly lower 
(higher) minority carrier diffusion lengths. 

At present, we can prove not even one of these four assumptions. In 



the following we will, nevertheless, discuss qualitatively what we expect for 
various experimental situations in terms of the model based on the above 
assumptions. 


4.2 Influence of a high temperature heat treatment 

A heat treatment at say 1200°C for 1 hour changes the microstructure of 
EFG ribbons. E.g., internal stresses in the ribbons may be minimized by 
movement of dislocations and their rearrangement into low-angle grain 
boundaries. We assume that, in addition, the self- interstitial content of 
grown-in carbon-self- interstitial agglomerates is decreased by dissociation of 
self-interstitials from the agglomerates and their subsequent diffusion to 
self-interstitial sinks, such as surfaces or dislocations. If we suppose that 
we start at the beginning of the heat treatment somewhere near the minimum 
electrical activity, then we expect to move to the left-hand side in Fig. 20 
and to get a decreased diffusion length as experimentally observed. 

A possible side-effect of such a high-temperature heat treatment is 
probably that dislocations which have been introduced by plastic deformation 
during the cooling process after crystal growth increase their ability to 
absorb intrinsic point defects as discussed in subsection 2.6. 


4.3 Interaction with metallic impurities ("Gleichmann Effect"). 

The carbon-rich zones produced during a high temperature heat treatment, 
as described in the previous subsection, are likely to act as gettering 



centers for metallic impurities. R. Gleichmann^ suggested that these 
submicroscopic agglomerates could possibly getter metallic impurities even 
more effectively than dislocation or high-angle grain boundary. This effect 
could explain why Gleichmann^ finds a much lower electrical activity 
(measured by EBIC) after a one hour 1200°C heat treatment than before the heat 
treatment: the metallic impurities which had decorated the dislocations 
before the heat treatment are redistributed to the carbon-rich zones during 
the heat treatment. It is an essential assumption in this context that the 
metallic impurities should, on the average, show a lower electrical activity 
when gettered in carbon-rich zones compared to the case when they decorate 
dislocations or grain boundaries. 


4.4 Influence of Phosphorus Didfusion 

Phosphorus diffusion starting from a high surface concentration, as 
generally used for producing the electrical p-n junction of a solar cell, 
generates a high supersaturation of silicon self- interstitials as discussed in 
subsection 2.4. Due to the high density of dislocations, which may act as 
sinks for self- interstitials, it can be expected that this high 
3upersaturation is restricted to some |jm near the phosphorus-diffused 
surface. The presence of a high supersaturation of intrinsic point defects in 
900°C phosphorus-diffused EFG ribbons can be also concluded form the electron 
microscopical observation of helical dislocations in the surface near 
region. 35 f 76 The 3U persaturation Cj/Cj 3 calculated from the number n^ of turns 
per unit length of such helices via eq (14) is in the order of 100-500, which 
can be expected for this kind of phosphorus diffusion at 900°C. 



Concerning the carbon-self- interstitial agglomerates, we assume that the 
high self-interstitial supersaturation will lead to an additional decoration 
of these agglomerates with self-interstitials. This means in Fig. 20 we move 
along the electrical activity curve from the left to the right. Suppose, we 
start near the minimum electrical activity, then during phosphorus diffusion 
the electrical activity will increase and the minority carrier diffusion 
length will decrease as often observed after a 900°C 30 min P-dif fusion. ^5 The 
situation is different if a high temperature heat treatment (e.g. 1 hr at 
1200° C) precedes the phosphorus diffusion, since then due to the high 
temperature treatment, we are likely to start to the left side of the minimum 
with a very small self interstitial content. In this case, the phosphorus 
diffusion will decrease the electrical activitiy (at least as long as the 
diffusion times are not excessive) and increase the minority carrier diffusion 
length as actually observed.? 5 

Phosphorus diffusion at higher temperatures than 900°C, say at 1050°C, 
produces a much lower self- interstitial supersaturation. Under these 
circumstances it can be expected that the interstitial content of the 
carbon-self-interstitial agglomerates cannot exceed a certain maximum value 
which should be lower than the corresponding maxiumum value for 900°C. 


5. Self-inter3titial induced oxygen precipitation in carbon-rich EFG ribbons 

It has been shown that phosphorus diffusion at 900°C for 30 minutes 
generates light-enhancement centers in EFG ribbons containing intermediate 
oxygen concentrations (several times 10^?cm”3) ,73»7^ suppose that these 
light enhancement centers are due to submicroscopic oxygen precipitates or 



basically the "new donors" discussed in- subsection 3 * 7 . Probably the 
nucleation of these Si 02 precipitates does not occur throughout the whole 
ribbons since the oxygen supersaturation is too low at 90 0°C. The nucleation 
obviously needs the triggering by phosphorus-induced silicon 
self-interstitials as also described in subsection 3 . 7 . Since this 
self-interstitial supersaturation is present only in a surface-near zone it 
can be expected that these light-enhancement centers are formed only in this 
surface-near region (whithin a few |jm of the surface) as is actually the 
case.' This is different in oxvgen-rich CZ crystals in which nucleation and 
growth of oxygen precipitates constituting the light enhancement centers may 
occur throughout the whole sample. The observation that no oxygen-related 
light-enhancement centers are formed any more at 1050°C ^5 is probably simply 
due to the fact that at this temperature the oxygen is no longer (or at least 
not sufficiently) supersaturated. 


6. Final remarks 

The speculative model presented for the influence of carbon on the 
minority carrier diffusion length in EFG ribbons certainly so far lacks 
experimental justification. The most important step would be to render 
visible the carbon-self-interstitial agglomerates. We suggest that copper 
decoration and subsequent X-ray topography would constitute a reasonable 
approach, in view of the fact that carbon- related B-swirls have also been 
detected this way.^ The self-interstitial supersaturaton generated by the 
diffusion of phosphorus also plays an important role in our model. If 
phosphorus as an n-type dopant were replaced by arsenic then we expect that no 
increase in the self-interstitial content of carbon-self-interstitial 



agglomerates would occur during diffusion (because As does not produce a 
supersaturation of self-interstitials) and all effects related to this process 
should disappear. Finally we mention again, that EFG ribbons are so 
inhomogeneous in their microstructure that we can well imagine that there may 
also exist regions, where the minority carrier diffusion length is not 
determined by carbon-self-interstitial agglomerates but rather by dislocations 
or metallic impurities as indicated in the "diffusion-length triangle" of Fig. 
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Figure 7 


Main factors likely to influence the minority carrier diffusion 
length Lp in EPG ribbons (diffusion length triangle) 

Point defects in silicon (schematically) 

Self-interstial (D^c®^) and vacancy (D^c®^) contribution to silicon 
self-diffusion according to G’cTsele and Tan]^ Experimental data 
points from references [ 16 — 19 3 - 

Crude estimate of vacancy (D y ) and self-interstial (Dj) 
diffusivities compared to other diffusivities in silicon? The 
uncertainty of the values of Dj and Dy is at least of the order of 
a factor of ten. 

Crude estimate of thermal equilibrium concentrations of vacancies 
(C«l) and of self-interstials (C® q )? the uncertainty of the 
individual values is at least of the order of a factor of ten. 

Self-interstial supersaturation = ( Cj - Cj^)/ C ® q as a function 
of KC1 content in an oxidizing atmosphere at 1200°C. Negative 
values of corresponde to a vacancy supersaturation] ^ 

a) Emitter-push effect 

b) braodening of a masked layer due to phosphorus-induced 



point-defect supersaturation (schematically). 
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8: Transmission electron microscope picture of a helical dislocation, 

which had formed below a phosphorus-diffused area. (Strunk et 
al 21 ) 

9: The concentration of A-type and B-type 3wirls in Ez silicon 

crystals versus growth rate (deKock, Rohsner and Booneu^jRohsner 
and Van den Boom^ 1 ). The physical relevant parameter is not the 
growth rate v itself, but rather the combination v/G, where 
G - I AT/ Ax I is the temperature gradient. 

10: Carbon in substitutional position is silicon(schematically ) 

11: Normalized change of lattice constant d of silicon due to 

substitutional carbon (C g ) or interstitial oxygen (C^) (according 
to references 42 and 50). 

12: Solid solubilities of substitutional carbon (C g ) an d of 

interstitial oxygen (0^) in silicon according to Beau and Newman^ 

and Craven 1 ^. 

13: Bond centered configuration of the oxygen interstitial in the 

silicon lattice. The numbers 1-6 note equivalent positions of 
the oxygen atom. A possible diffusion jump of the oxygen 

interstitial is indicated (according to Kaiser 1 ^, Hass 1 *? and 
Corbett et al^). 



Figure 14: Diffusivity of substitutional carbon (Cg)53 and of interstitial 

oxygen (0^) ^9*51 >52 j. n silicon as a function of reciprical 
temperature, compared to silicon self-diffusion (S^) , phosphorus 
diffusion and to the diffusion of interstitial copper (C.^). 

Figure 15: Configuration of the highly moblile carbon-self-interstitial 

complex (migration energy 0.9ev) found by Waters and Brower 51 * 
after low temperature electron irradiation. 

Figure 16: Volume condierations concerning the formation of SiC> 2 

precipitates, SiC predipitates or carbon agglomerates. 

Figure 17: Formation of stacking faults in carbon-rich silicon by phosphorus 

induced self-interstitial supersaturation ( Abe et al 50 ) 

Figure 18: Effect of increased nucleation of SiC> 2 precipitation by a 
phosphorus-induced supersaturation of self-interstitials 

(schematically, de Kock 5 ?). The carbon consectration of the 
sample is not known. 

Figure 20: Suggested electrical activity of carbon-self-interstitial 

agglomerates as a function of the self-interstitial content. The 
minimun electrical activity is expected for a stress-free 
configuration. It is also indicated by arrows in which direction 
the self-interstitial content shifts during a high temperature 
heat treatment or during phosphorus diffusion. 
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